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Abstract 

We  consider  solutions  of  the  Boltzmann  equation,  in  a  d-dimensional  torus,  <1  =  2, 3, 
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for  macroscopic  times  r  =  t/eN,  e  •«  1,  /  £  0,  when  the  space  variations  are  on  a  macroscopic  scale 
x  «=  eN~lr,  iV  j  2,  x  in  the  unit  torus.  Let  u(x,  t)  be,  for  t  g  t0,  a  smooth  solution  of  the 
incompressible  Navier  Stokes  equations  (INS)  for  N  =  2  and  of  the  Incompressible  Euler  equation 
(IE)  for  N  >  2. 

We  prove  that  (*)  has  solutions  for  t  g  t0  which  are  close,  to  0(e2)  in  a  suitable  norm,  to  the 
local  Maxwellian  (p/(2ir7’)‘,/J]exp(  -[y  -  eu(x,  t)]z/2T}  with  constant  density  p  and  temperature 
T.  This  is  a  particular  case,  defined  by  the  choice  of  initial  values  of  the  macroscopic  variables,  of  a 
class  of  such  solutions  in  which  the  macroscopic  variables  satisfy  more  general  hydrodynamical 
equations.  For  W  >  3  these  equations  correspond  to  variable  density  IE  while  for  N  -  2  they  involve 
higher-order  derivatives  of  the  density. 


1.  Introduction 


The  Boltzmann  Equation  (BE)  describes  the  time  evolution  of  the  positional 
and  velocity  density  of  a  gas.  An  old  question  is  the  relation  of  the 

Boltzmann  equation  to  the  Euler  and  Navier-Stokes  equations  of  fluid  dynamics. 
(For  an  historical  background  see  Chapman  and  Cowling  [1].)  In  particular, 
under  what  conditions  does  the  BE  lead  to  closed  hydrodynamical  type  equations 
for  the  locally  conserved  quantities  (collision  invariants)  which  depend  only  on  r 
and  r.  These  are  the  five  moments  of  /  corresponding  to  mass,  momentum  and 
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energy  densities, 


.  \ 

^  ' 
■ 


(1.1) 

s  -r 


p(r,  r)  =  f  f(r,  v,t)  dv,  pu(r,r)  =  f  v}(r,v,r)dv, 
pT(r,r)  =  \jd{v~  u(r,r))2f{r,v,r)  dv. 


The  possibility  ef,such  a  “contracted  description”  is  due  to  the  big  disparity 
(in  many  situatiohs)\>etween  the  microscopic  scale  of  the  BE,  represented  by  the 
mean  free  path  and  mean  free  time  and  the  macroscopic  scales  on  which  the 
(conserved)  hyditdyBamical  variables  change.  This  can  be  expected  to  lead  to  a 
state  close  to  locSal  equilibrium,  i.e.,  f(r,v,  r)  should  be  close  to  a  local 
Maxwellian, 


(1.2)  M(p,  u,  T\  v)  =  p/(2irT)d/2  exp{  -(v  -  u)2/2T), 


with  the  hydrodynamical  variables  changing  “slowly”  in  space  and  time.  M  is 
singled  out  by  the  fact  that  it  is  the  only  function  for  which  Q(M,  M )  =  0  (see 
[1H3]). 

The  disparity  between  microscopic  and  macroscopic  scales  is  the  central 
ingredient  in  the  derivation  of  a  variety  of  macroscopic  laws  from  microscopic 
model  systems  (see  [4]).  In  all  these  cases  one  follows  the  path  initiated  by 
Hilbert  [5]  and  Enskog  [6]  who  introduced  the  ratio  of  microscopic  to  macro¬ 
scopic  length  scales,  the  Knudsen  number  e,  as  an  explicit  parameter  in  the  BE 
changing  (*)  to 

(1-3)  dT/t(r,  t;  v)  +  v  Vr/t  =  fj. 


where  we  have  changed  the  order  of  the  variables  in  ft  to  emphasize  that  v  is  to 
be  thought  of  as  a  parameter. 

It  is  now  possible  to  make  various  formal  expansions  of  ft  in  powers  of  e 
about  some  local  or  global  Maxwellian.  In  particular  the  well-known  Hilbert- 
Chapman-Enskog  method  [1]  leads  to  a  sequence  of  equations  for  the  hydrody¬ 
namical  variables  corresponding,  respectively,  to  the  (non-dissipative)  Euler 
equations,  the  Navier-Stokes  equations,  etc.  From  now  on  we  always  consider  a 
collision  operator  Q  corresponding,  a  la  Grad  [2],  to  cut-off  hard  potentials. 

'A»-important  step  in  the  mathematical  development  of  the  theory  was  taken 
by  Caflisch  [7]  (and  by  Caflisch  and  Papanicolau  [8]  for  a  discrete  velocities 
model)  who  proved  the  following:  let  p,  u,  T  be  smooth  solutions  of  Euler 
equations  in  a  torus  for  t  e  [0,  /0J.  Then  there  exists  a  solution  fe(r,  r,  v)  of  the 
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BE  (1.3)  such  that,  for  t  <  t0,  we  have,  in  a  suitable  norm, 

(1.4)  \\ft-  M(p,u,T;v)\\<Ce. 

This  result  has  been  extended  by  Lachowicz  [9]  to  more  general  initial  data 
fe(r,  0;  v)  and  more  general  domains,  see  also  Bardos  [10]  and  Bardos  and  Golse 
[11]- 

We  now  observe  that  equation  (1.3)  is  a  special  case  of  scaled  equations  which 
can  be  obtained  directly  from  (*)  by  looking  explicitly  at  solutions  which  vary  on 
different  macroscopic  space-time  scales.  More  precisely  we  let  e-1  be  the  macro¬ 
scopic  scale,  i.e.,  the  size  of  the  domain  (torus)  in  which  the  system  is  confined 
and  we  let  r  =  e~lx,  t  =  e~at.  The  time  evolution  of  F(a)(x,  t;  v)  - 
/(e_1 x,  v,  e  at)  will  satisfy  the  (BE)0  equation 

dF‘ 

(1.5)  •  VxF('a)  =  F('J,  . 

with  some  initial  condition  F(a)( jc,0;  v).  Equation  (1.5)  is  nothing  more  than  a 
rewriting  of  (*)  in  terms  of  the  new  variables  x  and  t.  Any  new  results  valid  for 
“e  small”  must  therefore  follow  form  (*)  and  the  assumptions,  if  any,  on  the 
initial  state  F‘a)(xy0;  u).  This  will  be  discussed  further  in  Section  5;  see  also  [12]. 

Equation  (1.3)  thus  corresponds  to  a  =  1,  i.e.,  space  and  time  are  scaled  in 
the  same  way.  That  this  is  the  right  choice  for  the  Euler  equations  is  clearly 
related  to  the  fact  that  these  equations  are  invariant  under  this  scaling.  The 
reduction  to  Euler  type  equations  under  uniform  space-time  scaling  can  also  be 
proven  for  some  simple  microscopic  models  (see  [13])  and  shown  to  occur,  under 
the  assumption  of  local  equilibrium,  for  real  particle  systems  evolving  according 
to  Hamiltonain  equations  of  motion  (see  [4]). 

A  natural  question  to  ask  then  is  whether  there  is  any  different  micro-macro 
scaling  of  space-time  which  produces  the  Navier-Stokes  equations.  An  inspection 
of  those  equations  shows  that  they  are  not  invariant  under  any  rescaling  of 
space-time.  This  implies  that  there  are  no  space-time  scalings  which  produce  the 
full  Navier-Stokes  equations  as  a  well-defined  limit.  They  have  to  be  thought  of 
instead  as  a  correction  to  the  Euler  equations.  On  the  other  hand,  if  one  also 
scales  the  macroscopic  velocity  u,  i.e.,  one  considers  u ’s  very  small  compared  to 
microscopic  velocities  (or  the  Mach  number),  then  there  is  a  scaling  which  leaves 
a  “subset”  of  the  Navier  Stokes  equations,  and  in  particular  the  INS  equations, 
invariant;  c.f.  [14],  Similarly  there  are  scalings  of  u  which  make  the  dissipative 
term  in  the  INS  vanish  but  leave  a  “subset”  of  the  Euler  equations,  including  the 
IE,  invariant;  see  [15],  This  suggests  a  derivation  of  these  equations  directly  from 
the  BE  with  proper  scaling  of  the  initial  u’s.  We  carry  out  this  program  for  the 
INS  and  IE  in  Sections  2, 3,4.  We  leave  for  Section  5  the  discussion  of  the  larger 
class  of  hydrodynamical  evolutions,  obtained  by  imposing  some  constraint  on  the 
d  +  2  variables  occurring  in  the  compressible  hydrodynamic  equations,  which 
can  be  obtained  from  the  BE  by  choosing  suitable  initial  conditions. 
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2.  Statement  of  Results 

Let  u(x,  i)  be  the  solution  of  the  INS  or  IE  in  a  unit  {/-dimensional  torus  II, 
i.e.,  u  satisfies  the  equations 

(2.1a)  pd,u  +  pu  •  Vxu  =  -  Vxp  +  i)vAu, 

(2.1b)  div  u  =  0, 

(2.1c)  u(x,  0)  =  u0(x), 

where  p  is  the  constant  mass  density  and  i)N  is  the  viscosity  coefficient  for  N  =  2 
(INS),  while  ifN  =  0  if  N  >  2  (IE).  The  vector  field  Vxp  which  corresponds 
physically  to  the  pressure  gradient  is  not  specified  a  priori  but  is  to  be  considered 
formally  as  an  unknown  function  also  determined  from  (2.1).  The  parameter  N 
will  be  defined  later.  It  is  known  (see  for  example  [16])  that  if  u0  is  sufficiently 
smooth  there  is  a  t0  depending  on  u0  such  that  equation  (2.1)  has  a  unique 
smooth  solution  ( u(x ,  /),  P(x<  0)  f°r  t  I  to-  1°  two  dimensions,  t0  can  be  taken 
arbitrarily  large. 

In  order  to  define  precisely  the  smoothness  we  require,  we  introduce  the 
Sobolev  space  H^fll)  of  order  s  with  scalar  product 

(2.2)  (fi,g).=  t  /  dxv'h(x)-vJg(x). 

j- o'" 

Furthermore,  we  denote  by  BJS  the  space  of  the  measurable  functions  on 
n  x  R',  such  that 

(2.3)  11%,  =  sup  (1  +  v2)j/2\h(-,o)l,  \h\2s  =  (h,h)„ 

veKd 


is  finite. 

We  wish  to  prove  now  that  any  smooth  velocity  and  pressure  field  satisfying 
(2.1)  can  in  fact  be  obtained  from  the  solution  to  the  (BE)„  equation  (1.5)  for 
suitable  initial  distribution  and  choices  of  the  scaling  a.  The  case  a  =  1  is 
covered  by  the  Caflisch  theorem,  equt (1.4).  We  shall  now  consider  the  case 
1  <  a  g  2,  and  observe  that  by  letting  t*~'  *  e  we  can  rewrite  equation  (1.5)  as 
follows: 


dF‘ 

(2.4)  +  t~lv  •  VXF&  =  e-NQ(Fl,  F„), 

with  N  =  a/(a  -  1),  so  that  N  ^  2.  We  assume  that  N  is  an  integer:  our  proofs 
presumably  work  also  for  any  rational  number  N,  but  we  did  not  treat  that  case. 
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We  prove  the  following 

Theorem  I.  Given  a  divergence-free  field  u0  e  Hs(Tl),  s  ^  4,  p  >  0,  T  >  0, 
N  ^  2,  there  are  fixed  t0>  0  and  £0  >  0  such  that  equation  (2.4)  has  a  solution 
F„(x,  t\  v)  e  Bj2  for  t  g  t0,  j  >  4,  e  <  e0,  satisfying 

(2.5)  \\Fn(x,  f,  v )  -  M(p,eu{x,t),f-,  n)||>2  <  Ce2, 

whenever  u(x,  t )  is  a  solution  to  equation_{2.\)  in  HS(H)  for  t  g  t0.  In  the  case 
N  =  2,  the  viscosity  coefficient  tj  =  i}(p,  T)  is  given  by  the  equation  (3.32)  below. 

Remarks,  (i)  The  bounds  in  the  above  statement  are  not  optimal,  but  they 
are  sufficient  for  d  =  2,3. 

(ii)  Note  that  u  in  equation  (2.5)  is  not  the  same  as  the  u  in  equation  (1.1), 

U  — >  EU. 

(iii)  We  actually  prove  that  there  is  a  solution  Ff(x,  t;  v)  to  equation  (2.4) 
of  the  following  form:  we  consider  for  concreteness  the  case  N  =  2  and  write 
F{  =  F\ 

(2.6)  F'( x.  f,  v)  =  M(pe,  ut,  Tt)  +  e2g(x,  v,  t)  +  e2he{x,  t,  v), 
where 

(2.7a)  pt  =  p  +  e2p(2)(x,  t ), 

(2.7b)  ut  *  eu{x,  t)  +  e2u<2)(x, /), 

(2.7c)  Tt=  F+  e2T(2\x,t), 

g(x,  v,  t)  is  a  function  orthogonal  to  the  space  generated  by  the  collision 
invariants  (see  (3.25))  and  all  the  functions  pt,  ut,  Tt,  g(x,  v,  t)  and  he(x,  t,  v) 
are  uniformly  (in  e)  bounded  in  BJ2. 

Finally,  u(x,  t)_ in  (2.7b)  is  the  solution  to  the  INS  and  the  term  proportional 
to  e2  in  P'Te,  i.e.,  7 p<2)(x,  /)  +  pT<2)(x,  t)  is  “the  pressure”  p(x,  t)  associated  to 
the  solution  u(x,  t )  of  equation  (2.1)  with  N  =  2. 

(iv)  Our  results  involve  choosing  an  initial  distribution  Ff/(x,  0;  v)  which  is 
of  the  “right  form”.  This  may  in  fact  not  be  necessary.  As  was  shown  by 
Lachowitz  [9]  for  the  case  a  =  1  considered  by  Caffisch  (and  by  Ellis  [17]  and 
Ellis  and  Pinsky  [12]  in  some  simple  linear  case),  choosing  a  “slightly  wrong” 
distribution  at  t  =  0  has  the  effect  of  producing  an  “initial  layer”  which  then  has 
to  be  “matched”  to  the  hydrodynamic  solution.  Presumably  this  is  also  the  case 
for  a  >  1,  at  least  within  a  certain  class  of  distributions.  We  discuss  this  and 
other  questions,  including  extensions  of  Theorem  1  to  the  case  in  which  the  gas  is 
subject  to  a  suitably  scaled  “smooth”  external  force,  in  Section  5.  We  also  discuss 
there  choices  of  initial  conditions  which  lead  to  variable  density  IE  and  to  some 
“higher-order”  Navier-Stokes  equations  in  which  the  pressure  is  constant  up  to 
order  e2. 
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The  idea  of  the  proof  of  (2.5)  is  similar  to  that  of  Caflisch  [7];  perform  a 
truncated  (Hilbert)  series  expansion  in  e  and  control  the  remainder  equation.  All 
the  estimates  of  the  terms  of  the  expansion  and  of  the  remainder  are  due  to  the 
good  control  of  the  appropriate  hydrodynamic  equations  in  the  interval  [0,  /0],  In 
the  next  section  we  develop  the  expansion,  while  Section  4  is  devoted  to  the 
control  of  the  remainder. 


3.  The  Expansion 

We  first  consider  the 

Case  N  =  2.  In  this  case  equation  (2.1)  becomes 
(3.1)  d,F'  +  e~h}‘  vFe  =  £~2Q(Fe,  F'). 

We  look  for  a  solution  of  (3.1)  of  the  form 

(3.2a)  F'  =  mt+  L  e"/„  +  e*fR, 

n  —  2 


where 


(3.2b)  mt(x,  t,  v )  =  M(p0(x,  t ),  m(x,  t ),  T0(x,  t );  v). 


for  some  functions  u,  p0,  T0  to  be  determined.  The  functions  /„  will  be  chosen  to 
coincide  with  those  obtained  from  the  Hilbert  expansion  while  the  remainder 
term  fR  will  satisfy  the  equation  below.  We  note  that  mt  is  such  that  Q(me,  mt) 
=  0  and  that  the  first  correction  to  mt  is  0(e2). 

We  write 


(3.3) 
where 

(3.4) 


mt  =  m0+  £  e\  +  ek+1rk+1, 

n  —  1 


l  dn 
^  =  J?m'\ 


«-0 


n  1, 


and  rk  +  1  is  a  remainder  of  order  k  +  1, 

1  _  (v,Vj-  T&j) 

“  T*  lft  ~ 


(3.5) 


u  •  v 

Vi  ~  mo  f  >  92  ~  2m°~ 


To2 


UjUj,  i,  j  =  d. 


In  (3.5)  and  below  we  sum  over  all  repeated  indices. 
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Substituting  (3.2a)  and  (3.3)  in  equation  (3.1)  and  equating  terms  of  the  same 
order  in  e  not  containing  fR  we  get 


(3.6a) 

(3.6b) 

(3.6c) 


e_1:  v  •  Vm0  =  0, 

e°:  d,m0  +  v  •  V<f>!  =  Lf2, 

e1:  +  v  •  v(/2  +  <p2)  =  Lfi  +  L(l)f2 


and,  for  2  ^  Ac  ^  6, 

ek-  ^(/*  +  <P*)  +  wV(/*+i  +  «P*+i) 

(3'7)  =  Lfk+2  +  Lwfk+l  +  I  Q(fn,fJ. 

n,  m^2 
n+m-k+2 

Moreover  we  have  the  remainder  equation 

(3.8)  d, fR  +  e-h,  •  V/R  =  e~%fR  +  TjR  +  t2Q(fR,  fR)  +  e3A. 

In  equations  (3.6),  (3.7)  and  (3.8)  we  have  used  the  following  definitions: 


(3.9) 

LJ  —  2Q(me,  /), 

(3.1Q) 

L/=  2Q(m0,  /), 

(3.11) 

L(1)/=  2Q(rv  /), 

(3.12) 

zr/=2ef  fa"- 

\  n-2 

(3.13) 


A  -  -  d,[(/7  +  <p7)  +  e(/g  +  rg)]  -  v  •  v(/g  +  rg)  +  L(1)/g 

+  Le*“9  I  GU./J- 


*£9  n.mj  2 

n  +  m  —  k 


Now  we  analyze  the  consequences  of  equations  (3.6)  and  put  d  =  3:  equation 
(3.6a)  implies, 


Vpo  =  o,  vr0  =  0. 


(3.14) 
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Equation  (3.6b)  implies 

(3.15)  P(d,m0  +  v  •  V«pt)  =  0, 

(3.16)  (1  -  P)(d,m0  +  v  •  V<f>!)  =  Lf2, 
where 

(3.17) 

Pf(v)  =  f  f(v')dv'  +  v(  v'f(v')  dv'  +  (v2  -  3T0)  f  (v'2  -  3T0)/(v')  do'. 
Equation  (3.15)  thus  becomes 

(3.18a)  f  (d,m0  +  v  •  V<Pi)  dv  =  0, 

JR} 

(3.18b)  f  v(d,m0  +  v  •  V<p1)  dv  =  0, 

JR} 

(3.18c)  f  (v2  -  3T0)(d,m0  +  v  ♦  VqPj)  dv  -  0. 

•'r1 

Using  (3.14)  and  the  relations 

(3.19)  f  m0dv  =  p0,  f  V'VjiriQ dv  =  p0T08tJ,  (  v,q>l  dv  =  p0u,, 

J  R3  •'R3  •/R3 

(3.18a)  and  (3.18c)  give 

(3.20)  d,p0  +  di  vp0u  =  0, 

(3.21 )  3  p0d,T0  +  2  div  T0p0u  =  0. 

Integrating  equations  (3.20)  and  (3.21)  on  the  volume  II  it  follows  by  (3.14)  that 
d,p0  =  0  and  d,T0  =  0.  This,  together  with  (3.14)  implies  that 

(3.22)  p0  =  p,  T0  =  f, 
i.e.,  w0  is  a  global  Maxwellian, 

(3  23)  m0(x,  t;  v)  =  pg0  =  M(p,0,f;  v). 
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Moreover,  by  (3.22)  and  (3.20)  we  must  have,  for  the  expansion  to  be  consistent, 

(3.24)  div  u  =  0  (incompressibility). 

Finally  (3.18b)  is  satisfied  by  the  properties  of  the  Gaussian  integrals. 

Equation  (3.16)  yields,  for  all  /  ^  t0,  including  the  initial  distribution  at 
t  =  0, 


(3.25) 


fi  =  L~\T  l&o{vivj  ~  kv%j))pdjUi 

+  g0[p(2)  +  P(2)«(2>  *  o  +  em(v2  -  3f )], 


where  p(k\  u(k\  e(k)  will  denote  here  and  below  the  components  of  fk  along  the 
collision  invariants.  Since  they  are  in  the  null  space  of  L,  they  are  not  fixed  by 
equation  (3.25)  and  will  be  fixed  later.  We  find  (see  for  example  [3]) 

(3.26)  L~1^g0{vivj  -  $0%))  =  -*(t>2)(t>,”,  -  W 
with  'I'  a  known  positive  function.  Therefore, 

h  =  -&*(»z)(vlvJ-lvl»lj)godJul 

(3.27)  1 

+  g0[p(2)  +  p(2>«<2>  •  v  +  e(2)(v 2  -  3T)]. 

The  stress  tensor  and  the  heat  current  vector  associated  to  fk  are  defined  by 

(3.28)  ir<*>  =  /  vtOjfk  dv,  q]k)  -  /  v,(v2  -  3 f)fk  dv. 

J R3  ■/R3 

To  compute  the  ones  associated  to  f2  we  note  that 

/  <A;¥(p2)(p(i>  -  W^j)vkv,g0(v) 

(3.29)  R 

=  |i(f )(Mjf  +  «,A)(l  -  su)  +  ®( f)M«- 

with  positive  p(T)  and  o(T).  Moreover, 

/  P<2)uj\vkv,g0(v)  dv  =  0,  f  pmvkv,g0(v)  dv  =  p(2)f8kl, 

JR}  'R3 

j  emvkv,(v2  -  3 T)g0{v)  dv  =  2emT28k,  =  pTm8k,. 

•>  D3 


(3.30) 
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Hence,  using  (3.24), 

(3.31)  tt,(2)  =  -  \t](d,Uj  +  djU ,)  +  pSiJt 

where  the  viscosity  coefficient  and  the  correction  to  the  pressure  are  given  by 


(3.32)  ,  p  =  (p(2)  +  2Tem)T  =  p™T  +  pTm- 


On  the  other  hand,  since 

(3.33) 

(  rfo’^(u2)(o/o,  -  W&ij)viP2g o(v)  =  /  dvvkv2g0(v)  =  0, 

•'R3  •'R3 

and 

(3.34)  /  (u2  -  3T)v,Vjg0(v)  dv  =  2T28ip 

•'R3 

we  have 

(3.35)  qp  =  2r2p(2)u,(2). 

Using  the  relations 

(3.36)  f  <p2vdv  =  0,  f  <p2v(v2  -  3T)  dv  =  0, 

Jn  3  Jtt3 

and 

(3.37)  J"<p2viVjdv  =  puiuj, 


we  have  by  (3.6),  to  order  e, 

\ 


(3.38) 


div  p(2)u<2)  =  0, 


(3.39) 


d.pu  i  +  djirW  +  djpUjUj  =  0, 


h  =  L-‘[(  1  -  P)(3,<p1  +  u  •  v/2)  -  l<13/2] 


(3.40) 

+  g0[p(3)  +  P<3)«(3>  ’  o  +  em(v2  -  3 r)], 

with  p(3>,  m<3),  e<3)  to  be  fixed  later.  Using  the  expression  (3.31),  equation  (3.39) 
becomes  the  incompressible  Navier-Stokes  equation. 
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By  (3.7),  we  have,  for  2  ^  k  ^  6, 

(3-41)  p(d,(fk  +  <Pk )  +  v  V(A+1  +  <p*+i))  = 

fk  +  i  =  L  ^(l  —  p)[dt(fk  +  9t)  "*■  v  *  WA  +  i  +  <P/t  +  i)] 


(3.42) 


~L(l)fk+i  -  I  G(/.,/J 


n*  w ^2 


+  g0[p</£+2)  +  p(*+2)«(*+2)  •  y  +  et-k+2)(v2  —  371)]. 

After  equation  (3.41)  is  satisfied,  equation  (3.42)  determines  fk+2  up  to  the 
projection  on  the  space  of  the  collision  invariants.  On  the  other  hand,  we  can  use 
the  arbitrary  functions  p{k\  u(k\  e(k)  to  satisfy  equation  (3.41).  The  analysis  of 
the  conditions  for  p(k\  u(k\  e{k)  is  simple  and  we  omit  it:  the  functions  fk  can  be 
determined  and  will  be  considered  as  known  functions  of  x,  v,  t  below.  It  is  a 
classical  result  of  Grad  [2]  that  they  have  the  same  spatial  regularity  as  u  and 
satisfy  the  bounds 

(3.43)  I fk(x,  t,  u)|  £  9 k(x ,  t;  v)M0(v), 

where  the  & k(x ,  t;  v)  are  some  polynomials  in  v  with  coefficients  depending  on 
x,  through  u(x,  t). 

In  consequence  of  the  expansion,  if  u  is  a  solution  of  the  incompressible 
Navier-Stokes  equation  and  the  fk  are  chosen  as  above,  equation  (3.1)  becomes 
equivalent  to  equation  (3.8).  This  equation  is  much  simpler  than  the  original  one, 
since  the  nonlinear  term  is  multiplied  by  a  power  of  e.  This  fact  allows  then  for 
control  of  the  remainder  term. 

Case  N  >  2.  We  use  the  same  procedure  as  before.  We  look  for  a  solution  of 
(2.1)  of  the  form 


(3.44) 


2N  +  4 


Fl  =  mt+  £  en+Nfn+N+eN+2fR, 


where 


(3.45) 


me(x,  v,t)  =  M(pt(x,  t),  ut(x,  t),  T'(x,  t);v). 
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We  assume  that  pe(x,  t)  and  Tc(x,  t)  have  the  following  expansions: 

Pf(x,t)  =  p  +  £  pn(x,t)e", 

n- 2 

(3.46) 

\-i 

Tt(x,  t )  =  T+  £  Tn(x,t)e\ 

n  =  2 

and  «t  =  eu  for  N  =  3  having  the  expansion 

N-  1 

(3.47)  ue(x,t)  =  EU+  Y,un(x,t)En 

n- 3 


for  N  >  3. 

We  could  have  chosen  me  as  in  (3.45)  also  in  the  previous  case:  this  would 
have  only  changed  the  definitions  of  the  fk  in  (3.2)  for  k  ^  2. 

We  have  (see  equations  (3.4)  and  (3.5)) 

N  +  S 

(3.48)  m,  =  ffl0+  e<p!  +  £  En(<pn  +  Gn )  +  eN+6rN+6, 

n- 2 

where  we  denote  by  <p„  +  G„  the  coefficients  of  the  power  series  of  me,  while  the 
<pn  still  denote  the  quantities  (3.5).  We  now  proceed  as  in  the  case  N  -  2  and 
consider  the  equations  coining  from  the  orders  £-1,  e°  and  e1.  The  e- ‘-order  is 
exactly  the  same,  and  therefore  from  it  we  get  the  same  informations  as  before. 
The  e°-order  gives 

(3.49)  dtm0  +  v  •  =  Lf„. 

Equation  (3.49)  is  equal  to  (3.6b)  and  so  as  before  it  gives  the  incompressibility 
condition  and  the  expression  for  fN  as  in  (3.25).  The  e1- order  is 

(3.50)  dl9l  +  W(G2  +  Tl)  =  LfN+l  +  L»%. 

If  we  multiply  by  Vj  both  sides  of  (3.50)  and  take  the  {^-integral,  we  get  (see 
equation  (3.37)), 

(3.51)  pd,u  +  pu  •  Vxu  +  Vxp  =  0,  p  =  Tp2  +  pT2. 

Equation  (3.51)  is  the  incompressible  Euler  equation.  The  functions  /„  are 
determined  as  before  and  the  functions  p„,  Tn,  un,  n  =  3,-  •  •,  N  ~  1,  are  deter¬ 
mined  in  such  a  way  that  the  solvability  conditions  for  fN+n_x  are  satisfied, 


(3.52) 


+  G„-i)  +  o  •  V(<p,  +  G,)]  =0,  n  =  -  1. 
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Once  all  the  functions  appearing  in  the  expansion  are  determined,  we  are  left 
with  the  following  equation  for  the  remainder  fR. 

Let 


(3.53) 

and 


.  '2N  +  4 

LJ=2Q\  E  en~2fnJ 

h  —  2 


a  =  -d, 


2JV  +  4 

*Pa  +  S  ^iV  +  i  "h  £rN  +  6  "*■  E  £  fn 

n  *  N  +  5 


(3.54) 


-  v  •  V 


2N  +  4 

_  v  cn-N-(sf 

rN  +  6*  2-  6  Jn 

«  =  V  +  6 


+  £'%,* 4+  £  £  «/../.)• 


fcfi2V+5 


n, 


Then  /*  must  satisfy 

(3.55)  +  e-'u  •  V,/«  =  £-%/«  +  TjK  +  e2<2(/*,  /*)  +  e3,4, 

which  is  completely  analogous  to  equation  (3.8). 

4.  The  Remainder  Equation 

The  study  of  the  remainder  equation  follows  closely  the  works  [7],  [9].  There¬ 
fore  we  only  sketch  the  proofs  in  the  case  N  =  2.  The  main  difficulty  in  the  study 
of  equation  (3.8)  comes  from  the  negative  powers  of  e.  The  e_I  in  front  of 
v  •  vfR  may  be  easily  managed,  because  the  free  stream  operator  generates  an 
isometry  in  all  the  spaces  we  consider  below,  independently  of  e.  A  way  to 
control  the  e~2  in  front  of  Lt  is  to  use  the  nonpositivity  of  Lt  on 
L2[m;\x,  t;  v)  dv).  Following  Grad  [2],  one  would  consider  the  equation  for 
fR=  (mt)~l/2fR,  Since  mt  depends  on  jc  and  t,  the  equation  for  f R  would  contain 
polynomials  in  velocity  which  are  unbounded  for  large  velocities.  This  difficulty 
has  been  solved  by  Caflisch  [7],  by  means  of  a  decomposition  of  fR  in  low  and 
high  velocity  parts.  The  low  velocity  part  is  then  controlled  by  the  usual  Grad 
technique,  while,  to  control  high  velocities,  one_introduces  a  global  Maxwellian 
M*(v)  with  a  temperature  T*  larger  than  T,  in  order  to  gain  from  it  a 
convergence  factor  in  the  velocity  space. 
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More  precisely,  following  Caflisch  [7],  we  look  for  the  solution  of  equation 
(3.8)  in  the  form 

(4.1)  /*  =  fag  +  jATh, 
where 

(4.2)  M*{u)  =  M(p,0,T*,v),  T*  >  f, 

and  g  and  h  (the  low  and  high  velocity  parts,  respectively)  are  defined  as  the 
solutions  of  the  following  coupled  system  of  equations: 

(4.3)  (d,  +  e~1v  •  V  )g  =  e~2&tg  +  e^xo'^'h, 

(d,  +  £~xv  -v)h=  -  pog  -  £~2{ve  -  x#*)h  +  Le{ag  +  h) 

(4.4) 

+  eVr*(ag  +  h,  ag  +  h)  +  e3a, 

(4.5)  g(x,  u,0)  =  h{x,  u,0)  =  0, 
where,  given  y  >  0, 


lh 

M  ^  y. 

(4.6) 

x(«)  = 

X  =  1  -  X- 

(o 

otherwise. 

(4.7) 

M  =  \me  '{<>1  +  e  lv  v)me. 

a  ~  M*  ’ 

(4.8) 

J* 

II 

> 

dn(vl  -  v) 

•  n  mt( x,  t,  w, 

•'RJ  *'(ol  —  v) 

In^O 

(4.9) 

v*  =  f  dvl  f 

dn(vl  - 

v)  •  n  M*(vx) 

(4.10) 

^e/=(K) 

’i.ift/i  -  - 

-V'f+JTJ, 

(4.11) 

(4.12)  J?/«  (yfKF)~l  Lt[y[Kr  f\, 

(4.13)  P*r*(/,g)  =  {yfW)~lQ[y[XFf,y/Wg\, 

A 


(4.14) 
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The  term  n  in  equation  (4.4)  is  unbounded  in  the  velocity,  but  now  it  is 
controlled  by  the  exponentially  decaying  factor  a. 

To  study  this  problem  we  consider  a  linear  version  of  the  equations  (4.3)-(4.5), 
i.e.,  we  replace  g)  +  ea  by  a  given  function  b.  This  yields 


(4.15)  (d,  +  e~lv  •  v)g  =  e2£Ttg  +  £  2X° 


(4.16) 


( d,  +  e  •  V  )  ft 

=  —f nag  -  £-2(vt  -  +  J?{og  +  h)  +  e2b. 


(4.17) 


g(x,  v , 0)  =  h(xy  u,0)  =  0, 


Now  we  have  a  system  quite  similar  to  (6.1)-(6.3)  of  [7].  The  main  difference  is 
that  the  factor  e_1  in  the  right-hand  sides  of  (6.1)  and  (6.3)  of  [7]  is  replaced  by 
e~2  and  that  the  left-hand  sides  of  (4.15)  and  (4.16)  contain  the  factor  e-1. 
Moreover,  the  x  variable  is  three-dimensional  instead  of  one-dimensional.  This 
requires  the  control  of  higher  derivatives  in  x,  because  eventually  the  nonlinear 
term  has  to  be  bounded  using  the  algebraic  properties  of  Jf?s  fors>  \d.  The  e~2 
factor  in  front  of  is  not  relevant  because  the  proof  is  based  on  the 
nonpositivity  of  .£?.  Finally  the  factor  e_1  in  front  of  v  *  v  does  not  make  any 
problem  as  already  mentioned  (see  Appendix).  We  therefore  prove  the  following 

Lemma.  There  is  an  e0  >  0  such  that  any  solution  to  equations  (4.15),  (4.16) 
and  (4.17)  satisfies  the  bounds 


(4.18) 

(4.19) 


llffll js  S  e1/4||  b/V'\\j+2„ 

PIU  ^  e5/4PA IU. 


for  j  >  3  and  s  ^  2  if  e  <  e0 


The  proof  of  the  lemma  is  almost  a  repetition  of  the  proof  in  [7],  [9];  it  is 
given  in  the  appendix  for  the  sake  of  completeness. 

Now  we  assume 


(4.20)  b  -  v*T*(og  +  h,  og  +  h)  +  ea. 

Since  in  three  dimensions  the  space  ^  is  an  algebra  for  s  >  |,  using  the  Grad 
estimate  for  T*, 


(4.21) 


lir*(/,  g)|| j3  <>  c\\f\\js\\g\\j,  for  j  >  I,  s  >  §, 
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and  the  bound  (3.43)  to  estimate  a,  we  get,  for  s  >  \,  j  >  \, 

(4-22)  | b/vXs  ^  (llfflloj  +  11*11?,  +  «). 

where  we  have  taken  into  account  the  exponentially  decaying  factor  a.  Therefore, 
by  the  lemma, 

(4-23)  l|*l|>+2,  g  e5/4(llg|lL  +  ll*lly+2,  +  «) 

which  implies 

(4.24)  \\h\\j+2,  s  llg||0j  + 


for  c  small  enough.  Then, 

(4.25)  ||g||,,  <;  £1/4(||*||L  +  \\h\\2J+2s  +  ce)  ^  3^{\\g\\2s  +  cs) 

and  hence  the  estimate  for  *.  Using  this  and  (4.24)  we  get  the  estimate  for  h, 
(4-26)  11*10,  ^  Ce,  11*11,,  <;  Ce, 

for  some  constant  C. 

A  standard  iteration  procedure  gives  the  existence  of  the  solution  once  the 
a  priori  estimate  (4.26)  has  been  established.  This  proves  the  theorem. 

5.  Generalizations  and  Discussion 

Our  main  result.  Theorem  1,  can  be  restated,  see  comment  (iii)  in  Section  2,  as 
follows:  Let  F„(x,  r,  v)  be  the  solution  of  the  BE  (2.4)  for  the  initial  distribution 

2N  +  4 

(5.1)  /^(a,0;  v)  =  M(p,  eu0(x),  T;  v)  +  £  enhn(x,v), 

n  =  N 

(5.2)  divuo  =  0, 

where  the  h„(x,  v),  pt(x,0),  ut(x,0),Tt(x,0)  are  determined  from  the  initial 
values  of  the  hydrodynamical  variables  satisfying  (2.1).  (E.g.,  for  N  =  2,  h2(x,  v ) 
is  given  by  the  right  side  of  (3.27)  for  t  =  0  and  pe,  ue,  Tt  are  given  by  initial 
values  of  (2.7)  for  N  =  2  and  by  (3.46), (3.47)  for  N  >  2.)  Then  F„(x,  t ;  v) 
satisfies  (2.5)  for  t  <;  t0.  General  theory  assures  (or  ought  to)  that  this  solution  is 
unique. 
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A  natural  question  then  is,  what  happens  if  we  replace  (5.1)  by 
F*(x,0;  v)  =  M(p0(x),  eu0(x),  T0(x);  v ) 

(5.3)  2N+A 

+  eh1(x,v)+  X!  e"h„(x,v), 

n-2 

and  drop  (5.2b)?  The  answer  turns  out  to  be  rather  complicated  and  more 
dependent  on  N  in  (2.4)  than  before.  Let  us  first  state  the  result  for  the  case 
N  4.  We  assume,  instead  of  (3.44),  a  solution  of  the  form 

2N+4 

(5.4)  FZ-m,+  £  e”/„  +  tN+2fR. 

n- 1 

In  view  of  (3.46)  and  (3.47)  there  is  no  loss  of  generality  in  assuming  that  the  /„ 
are  orthogonal  to  the  collision  invariants  for  n  =  1,-  •  •,  N  —  1.  The  method  of 
Section  3  then  allows  us  to  prove  that  Ff,  is  a  solution  to  the  (2.4)  if  /„  =  0  for 
n  =  1,*  •  •,  N  -  2,  and  fN_l  satisfies  the  equation 

(5.5)  v  •  vm0  =  LfN_v 

The  compatibility  condition  for  (5.5)  implies  that  p0  and  T0  may  depend  on  x, 
but  the  pressure  p0  =  p0T0  has  to  be  spatially  constant.  Therefore  one  can  prove 
that  Theorem  1  is  still  true,  but  (2.5)  is  replaced  by 

(5.6)  ||F^  -  M(p0(x,  t),  eu(x,  /),  T0(x,  t );  i>)||  g  Ce2, 
with  p0  and  u  satisfying 


d,p0  +  u  •  Vp0  =  0, 


(5.7)  div«  =  0, 

p0(3,u  +  u  •  Vu)  =  -Vp, 

which  is  the  IE  with  nonconstant  density.  Moreover,  T0  is  given  by  p0T0  =  p0  for 
some  p0  >  0. 

It  is  clear  from  (5.5)  that  fN_1  depends  on  vp0-  In  particular,  for  N  =  3,  f2 
does  not  vanish  and  (3.5)  has  to  be  modified  accordingly.  However,  it  turns  out 
that  /2  is  an  odd  function  of  v  and  the  stress  tensor  7r(2)  vanishes.  Therefore,  the 
hydrodynamical  equations  are  the  IE  with  nonconstant  density,  also  in  this  case. 
Considering  now  the  case  N  =  2,  we  find  that  fx  does  not  vanish  and  (3.21), 
(3.27)  and  (3.39)  have  to  be  modified.  The  hydrodynamical  equations  which  can 
be  obtained  along  the  same  lines  as  in  Section  3  become  more  complicated  and 
involve  additional  terms,  containing  higher  derivatives  of  p0,  which  depend  on 
extra  transport  coefficients.  We  do  not  go  into  details  here  but  simply  write  them 
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down  and  then  make  some  remarks: 

(5.8a)  '  d,p0  +  div  p0u  =  0, 

(5.8b)  div(^p0r0«)  =  dt(*tJ{T0)djT0)t 

3,p0ui  +  dj[pAj  +  Po«,«y] 

(5.8c)  =  dj(n(T0)(d,Uj  +  dju,)  +  y(T0)8u  div  a) 

+  Bj{cLijkAT0)dkd,T0  +  Pijki(T0)ukd,Ta  +  Y,y*/d*7od/To)> 

(5-8d)  T0=g, 

with  k,  p,  p,  a,  P,  y  functions  of  T0  which  depend  on  the  intramolecular  interac¬ 
tions,  while  pQ  does  not  depend  on  x  or  t. 

Remarks,  (a)  If  the  density  and  temperature  of  the  initial  distribution  are 
constant  up  to  terms  of  order  e,  then  we  get  the  situation  discussed  in  Section  3. 
The  same  is  true  if  we  consider  isothermal  or  isoentropic  flows. 

(b)  The  scaling  considered  in  this  paper,  for  N  =  2,  leaves  the  new  equa¬ 
tions  invariant.  This  means  that  the  extra  terms  cannot  be  removed  by  scaling. 
On  the  other  hand,  if  one  assumes  that  k  vanishes,  then  (5.8b)  implies 


div  u  =  0. 

If  we  assume  that  v  does  not  depend  on  T0  and  that  a,  P,  and  y  vanish,  we  have 

div  u  =  0, 

d.Po  +  «  •  Vp0  =  0, 

p0d,u  +  p0u-  Vu  =  -  VPz  +  p&u, 

which  represent  the  INS  with  nonconstant  density. 

There  are  various  directions  in  which  the  present  paper  can  or  ought  to  be 
extended.  We  discuss  some  of  them  here. 

5.1.  Initial  and  asymptotic  distributions.  It  would  be  desirable  to  prove  a 
theorem  of  the  following  type.  Given  a  “smooth”  initial  distribution  f{r,  u,0), 
the  BE  (*)  will  bring  /,  in  times  short  compared  to  e-1,  to  a  local  equilibrium 
state  with  parameters  p,  u,  T  satisfying  the  compressible  Euler  equations  for 
times  e-1  <  r  <  e~“,  a  >  1,  the  IE  equations  for  times  r  -  e~°,  1  <  o  <  2,  and 
finally  the  INS  for  times  r  -  e~2  leading  eventually  as  t  -»  oo  to  the  global 
Maxwellian  distribution  in  which  p,  T,  and  u  are  constants  determined  by  the 
initial  state  (u  can  be  taken  to  be  zero  by  a  Galileian  transformation). 
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As  already  mentioned,  a  first  step  in  that  direction  was  taken  by  Lachowitz  [9] 
who  extended  the  Caflisch  theorem  to  the  case  when  the  initial  distribution  is  not 
“too  far”  from  a  local  Maxwellian.  His  result  for  the  BE  (1.5)  with  a  =  1  can  be 
stated  as  follows: 

Theorem  2  (Caflisch-Lachowitz).  Let 

HF^x.O;  «)  -  Af(p0,  u0,  T0;  i>)||  <  8, 

for  some  8  sufficiently  small.  Then  there  is  a  e0  >  0  s.t.  for  any  S'  >  0 

IIFJjCx,  t;  v)  -  M(p,  u,  T\  d)||  ^  Ce  for  all  e  <  e0, 

for  all  t  G  (S',  r0],  where  p(x,  t ),  «(x,  /)  and  T(x,  t)  are  the  smooth  solutions  of 
the  compressible  Euler  equations  up  to  time  t0  with  initial  values  p0(x),  u0(x),  T0(x). 

We  hope  that  the  work  in  progress  by  Bardos,  Golse  and  Levermore  [18] 
which  utilizes  the  recent  work  of  di  Pema  and  Lions  [19]  proving  the  global 
existence  of  weak  solutions  of  the  BE  will  in  fact  lead  to  results  along  the  lines 
described  above.  For  the  moment  however,  it  should  be  emphasized,  we  have  no 
uniform  bound  for  arbitrary  times  t  even  for  the  restricted  class  of  initial 
distributions  considered  here.  We  expect  however  and,  in  cases  where  the 
deviation  from  a  global  Maxwellian  is  small,  can  prove  that  as  t  -»  oo  the 
solution  of  the  INS  (2.1)  and  of  the  BE  (2.4)  are  u(x,  t)  -*  0  and  F\x,  f,  v)  -* 
Af(p,0,  f\  v)  as  t  -»  oo.  Here 

f  =  f  +  dvv2[F‘(x, 0;  v)  -  M(p,  eu0(x),  f;  i>)]  +  e2pul(x)  j 

is  determined  by  the  conservation  of  energy,  p  and  T  being  the  initial  density 
and  temperature. 

5.2.  Fluid  subject  to  external  forces.  The  systems  considered  so  far,  ap¬ 
proach,  as  /  -*  oo,  the  state  of  global  equilibrium  with  uniform  p,  u  and  T.  We 
wish  now  to  consider  the  more  interesting  situation  in  which  the  hydrodynamical 
flow  is  subject  to  an  external  force,  i.e.,  when  (2.1a)  is  replaced  by 

(5.9a)  pd,u  +  pu  •  Vu  =  -  Vp  +  qNhu  +  pb{x), 

where  b(x)  is  smooth  and  periodic.  (We  consider  for  simplicity  only  the  case 
N  =  2,  corresponding  to  the  INS.)  This  corresponds  to  changing  (2.4)  to 

<9  Fl 

(5.9b)  +  e~'v  •  V XF‘  +  eb  •  VvF'  =  e~2Q(F',  F‘). 

Note  that  in  the  original  Boltzmann  equation  (*)  the  force  term  would  be 
e3b(er)  •  vvf(r,v,  t),  i.e.,  the  force  varies  on  a  macroscopic  scale  and  is  very 
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small  compared  to  the  interparticle  forces  in  the  gas  which  give  rise  to  the 
collision  term.  It  is  now  straightforward  to  extend  Theorem  1  to  this  case. 
Equations  (3.2a)  to  (3.2b)  are  unchanged,  while  (3.6c),  (3.1)  and  (3.8)  are 
modified  in  obvious  ways.  The  only  difference  is  in  the  proof  of  the  lemma  (see 
Section  4)  given  in  the  appendix.  We  need  to  replace  the  free  motion  (x,  v)  -* 
(x  +  vt,  v)  by  the  evolution  in  the  force  field 

(x,  v)  -»  lx  +  vt  +  f'dsf  ds'b(x(s’),v  +  f'dsb(x(s)), 

\  Jo  Jo  Jo 

which  of  course  changes  also  the  various  constants  in  the  bounds. 

There  are  two  simple  special  cases  which  are  worth  considering: 

(a)  The  conservative  case.  We  take  b  =  —VxU(x).  In  this  case  the  energy 
of  a  gas  particle  \v2  +  e2U(x)  is  conserved  and  the  equilibrium  state  of  the  gas  is 
M(pe(x),0,  T;  v)  with  the  density  pe(x)  given  by 

(5.10)  pt(x)  =  pexp{-e2U(x)/2T  }. 

We  could  therefore  replace  p  by  pe(x)  in  (2.5)  and  state  the  theorem  in  terms  of 
this  local  Maxwellian. 

(b)  The  uniform  case.  We  take  b  =  i  independent  of  x  and  consider  an 
initial  distribution  F'(x,0;  v)  =  Af(p,  eit,  T\  v),  where  u  is  a  constant  (which  as 
already  noted  can  be  set  equal  to  zero  by  a  Galileian  transformation).  The 
solutions  of  the  INS  and  of  the  BE  (5.8)  and  (5.9)  are  then 

(5.11a)  w(x,  t)  =  u  +  bt, 

(5.11b)  F'(x,  /;  v)  =  M(p,  e[u  +  bt],  f;  v), 

which  are  seen  to  satisfy  (2.5)  with  C  =  0. 

The  interesting  case  is  of  course  when  b  *  -  vt/  but  fnbdx  =  0.  We  expect 
then  that  (5-8)  will  have  a  stationary  solution  «,(x)  which  will  be  approached 
asymptotically  as  t  -*  oo  starting  with  any  u„.  The  gas  distribution  F'(x,  t;  v)  on 
the  other  hand  does  not  approach  a  stationary  state  in  this  case.  The  reason  for 
this  is  that  the  external  force  continues  to  do  work  on  the  system  even  in  the 
stationary  state  which,  since  the  system  is  finite  for  e  #  0,  leads  to  an  ever 
increasing  temperature.  This  however  does  not  appear  at  the  level  of  the  INS 
(5.8a)  which  only  describes  the  organized  flow  field. 

To  see  this  explicitly  we  consider  the  case  where  u0(.x)  =  «,(x).  Multiplying 
(5.8a)  by  u,(x),  integrating  over  the  torus  II  using  (5.8b),  we  get 


(5.12) 
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On  the  other  hand,  if  we  multiply  (5.9)  by  \v2  and  integrate  over  both  v  and  x 
we  find 

(5.13) 

j  dx  J  dv  \v2F'(x,  f,  v)  =  p(v)2  +  \dpT]  =  £  JnP(v)  *  bdx > 

so  that  for  (v)  =  eus  the  rate  of  heat  production  in  the  system  is  proportional  to 
e2  times  the  right-hand  side  of  (5.12). 

It  appears  therefore  that  in  order  to  get  a  stationary  distribution  for  the 
positions  and  velocities  of  the  gas  particle  in  the  presence  of  a  stationary 
macroscopic  flow  one  needs  either  (i)  provide  contact  with  external  heat  reser¬ 
voirs,  e.g.  through  contact  with  walls  at  fixed  temperatures,  which  will  absorb  the 
energy  or  (ii)  consider  an  infinite  system  in  which  any  finite  heat  generation  can 
“dissipate  to  infinity”.  An  example  of  case  (i)  is  the  study  by  Bardos,  Caflisch 
and  Nicolaenko  [20]  of  the  stationary  distribution  of  a  gas  between  two  plates 
kept  at  different  fixed  temperatures.  An  example  of  the  case  (ii)  is  the  study  by 
Asano  and  Ukai  [21]  of  the  flow  past  an  obstacle.  They  prove  the  existence  of 
stationary  distributions  which  are  actually  attracting  for  a  certain  class  of  initial 
distributions. 

53  Linear  Boltzmann  and  incompressible  Stokes  equations.  Many  of  the 
questions  raised  in  the  last  sections  become  easier  to  answer  if  instead  of  starting 
with  (*)  one  starts  with  the  linearized  BE  for  /, 

(5.14)  d,f  +  v  •  Vxf  —  If , 
where 

(5.15)  Lf  =  2g[m0,  /], 

mQ  =  M(p,  0,  T;  v)  being  the  global  Maxwellian  with  density  p  and  temperature 
T.  Scaling  variables  like  in  Section  2  and  assuming  N  =  2  we  get,  instead  of  (2.4), 
the  equation 

are 

(5.16)  -jj-  +  VXF'  =  £~2LF'. 

In  this  case,  using  the  arguments  of  Section  3,  it  is  quite  easy  to  prove  the 
following 

Theorem  3.  Given  a  divergence-free  field  u0  e  tf,(II),  s  ^  2,  p  >  0,  T  >  0, 
equation  (5.16)  has  a  unique  solution  F'(x,  f;  u)  e  L2[mo\v)  dv  dx]  satisfying , 
for  any  finite  t. 


(5.17) 


||Fe(jc,  t\  v)  -  M(p,  eu(x,  t),  T;  i>)||  <  Ce2, 
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whenever  u(x,  t)  is  a  solution  in  J^(II)  of  the  Stokes  equation : 

(5.18a)  pd,u  =  -Vp  +  tjA  u, 

(5.18b)  div  u  =  0, 

(5.18c)  u(x,  0)  =  u0(x), 

and  ij  is  determined  as  in  Theorem  1. 

To  prove  Theorem  3  it  is  enough  to  realize  that  following  the  lines  of  Section 
2,  equations  (3.6)  are  modified  as  follows: 

(5.19a)  e~l:  v»vm0  =  0, 

(5.19b)  e°:  d,m0  +  v  •  VqPi  =  L(f2  +  <p2), 

(5.19c)  e1:  +  v  •  v(/2  +  <p2)  =  L{f2  +  <p3). 

From  (5.19b)  we  get  that  f2  =  f2  -  <p2,  where  f2  is  equal  to  the  right-hand  side  of 
(3.27).  Therefore,  (5.19c)  becomes 

(5.20)  d,<p2  +  v  •  V/2  =  L(/3  +  qp3). 

From  (5.20),  equations  (5.18)  easily  follow.  The  equation  for  the  remainder  is 
linear  and  this  allows  us  to  control  it  by  the  nonpositivity  of  the  operator  L  on 
L2lmo  l(v)  dvdx). 


Appendix 

We  use  the  following  estimates  proved  in  [2],  [7]  and  [9]: 

(A.l)  m\\JS  ±  c\\f\\j.u,  llJTVIly,  S  41/11, -u  for  j^hs^o, 

(A.2)  IWlloo  £  411/111,  llJTVIloo  *  411/111. 

0  <  V0  £  vt(x,  V,  t)  S  *l(/  +  \v\Y,0  <V0£  v*(v)  g  ^(1  +  |o|)*, 

0</JSl, 


(A.3) 

(A.4) 


£  c\\f\\Js, 

*  i > 


(A.5)  W\\  g  |^(v*JTe  )/|o  g  e‘||/||y+. 


BOLTZMANN  EQUATION 


1211 


In  (A.2)  the  |||  •  |||  denotes  L2-norm  in  velocities  and  positions.  Estimate  (A.4) 
follows  from  the  bounds  (3.43)  and  (A.5)  follow  from  a  direct  computation  using 
the  fact  that  the  temperature  is  constant. 

We  denote  by  x  +  vt  the  translation  on  the  torus  FI  with  velocity  v  in  time  t. 
Of  course  the  Lebesgue  measure  is  invariant  under  this  translation.  We  write 
(4.15)-(4.17)  in  integral  form: 


h(x,  v,  t)  =  jf'dsjexpj- J^ds'vt(s')e  2j 


•[-liag  +  e2xJf*h  +  Sf’iag  +  h)  +  B2b]\  , 


g(x,v,  0  =  jTds  jexpj  -  f‘ ds'vt(s')t 


•[e  2Jftg  +  e  2xo  , 


where  the  symbol  {  }#  means  that  the  functions  in  braces  have  to  be  computed  at 
(x  —  (t  —  s)ve~l,  v,  s).  Since  we  compute  integral  norms  in  x,  the  presence  of 
e_1  makes  no  difference  for  the  invariance  of  the  Lebesgue  measure  on  the  torus. 
Perform  the  time  integral  which  gives  a  factor  e2/vt.  Using  (A.l)  (A.4)  and 
choosing  y  large  enough  we  get,  by  (A.6), 

(A.8)  \\h\\j0  g  e2||g||00  +  e4  jbl  . 

'  llyo 

Differentiating  (4.16)  and  writing  it  in  integral  form  we  have 
Vh(x,  v,  t)  =  fds  jexp|  ~  J' ^s'vXs')e~2^j 
(A.9)  *[-v(n<7g)  +  e~2(v»t)h  +  e~2xV(  Jf*h) 

+  v(^(ag  +  A))  +  r2vf>]J  . 

Using  (A.8)  and  (A.l),  (A.2),  (A.4),  (A.5)  we  get,  for  y  large  enough, 


(A.10) 


IIV/ilU  *  ^llsllot  +  e4  Vb  ■ 
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By  the  same  argument  higher  derivatives  may  be  estimated  and  we  get 

(A. 11)  11*11,-,^  e2|l*llo,  +  e4|^|  • 

Now  we  estimate  g.  Following  [2]  we  take  advantage  of  the  bound  (A.l)  to 
reduce  the  estimate  of  ||g||>s  to  an  L2  estimate  which  follows  by  the  nonpositiv¬ 
ity  of  the  operator  Namely,  denoting  by  C(y)  a  constant  diverging  with  y,  we 
have 

1lgll,0  ^  m\jo  +  ^  Il^«lly0  +  C{y)\\X*h\\j0 

(A  12)-  n«ll>-i°  +  C(y)P||7_io  ^  llg||,-io  +  c(v)|£2llglloo  +  E“|  iqJ 

^  iigiiy-io + • 

II  e  II/- 10 

for  e  small  enough.  In  the  same  way,  using  (A.2)  we  get 

(A. 13)  llglloo  ^  111*111  +  £?/2  ■ 

*  00 

A  recursive  argument  (Grad  hierarchy)  then  gives 

(A.14)  l|g||,-03  |||g|||  +  e7/2)l yb\  . 

II  »  II/- 10 

The  L2  estimate  is  obtained  using  integration  by  part  to  cancel  the  streaming 
operator  e  ~  lv  •  V,  and  the  nonpositivity  of  i£t  to  get 

d,|||*|||  ^  £-2|||x®“1^III  ^  C(y)e_2|||jr*^j||  ^  C(y)ir2|||/.|||  <;  C(y)e-2P||20 

(a.is)  r  1  II  1  IM  1 

^  C(y)  Hglloo  +  *2  ^  C(y)  |||g|||  +  e2  -h  . 

L  «  II20J  II  e  20. 

Therefore,  by  the  Gronwall  inequality  we  have 

Hlglll  £  exp{C(y)/0}e2||^h||  ^  e3/2|^-6 

II  »  II20  II  * 


(A.16) 


20 
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if  e  is  small  enough.  Therefore, 
(A.17)  ||g||,0  ^  e3/2 


11,-10 


j>  3. 


The  equation  for  the  derivatives  of  g  can  be  studied  in  the  same  way.  Higher 
norms  of  b  appear  in  consequence  of  (A.5). 
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